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Abstract 

A theory of nonlinear convective depletion is set up as a nanosphere translates fast through 
a semidilute polymer solution. For nanospheres a self-consistent field theory in the Rouse 
approximation is often legitimate. A self-similar solution of the convective depletion equation 
is argued to be feasible at high velocities. The nature of the thin boundary layer in front of 
the propagating particle is analyzed. One example of convective depletion is when a charged 
protein moves through a semidilute polymer under the influence of a high electric field. The 
protein velocity is then proportional to the fifth power of the field. The theory could be useful in 
interpreting the separation of protein mixtures by microchip electrophoresis. 
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1 Introduction 



The propagation of colloidal particles through polymer solutions presents us with math- 
ematical problems of considerable complexity, both at low and high velocities. A major 
simplification results if we focus on a particle of nanometer dimensions for the following 
reasons: 

1) If the size is much smaller than the correlation length of, say, a semidilute polymer 
solution, the equilibrium depletion theory simplifies considerably 

2) If the solvent is intermediate rather than excellent as is often the case in practice (0), 
it suffices to use a self- consistent field description based on chain entropy alone The 
polymer may be viewed as quasi-ideal. 

3) It can be argued that the motion may be split up into two types of dynamics: dynamic 
depletion as a result of the rearrangement of the local polymer density about the particle, 



and the translation of the dressed particle (i.e. the partic 



^together with its depletion layer) 
). This is valid because the two 







through the strongly fluctuating polymeric suspension 
respective time scales are well separated. 

4) The polymer stress is highest close to the moving nanoparticle (jsl) so it is reasonable 
to suppose the dynamics of a chain section nearby is Rouselike or freely draining 

In effect, the polymer density is low in the depletion layer, for the density must be zero at 
the particle surface. 

5) In general, at very high velocities, the elasticity of the polymer must come into play 
which could lead to complex phenomena like shock waves An important dimensionless 
variable is then the ratio of the particle velocity to that of the elastic waves. But it is 
probably safe to neglect these when the size of the nanoparticle is much smaller than the 
distance between entanglements in the polymer solution. 

Besides being of intrinsic interest from a theoretical point of view, the fast translation of 
small particles through semidilute polymer does occur in practice, namely during capillary 
gel electrophoresis of biomacromolecules at high electric flelds ( jsl ) • 

I focus on the fast propagation of a nanoparticle, a problem reminiscent of others in 

n 



condensed matter physics (e.g. convective diffusion ( runaway electrons in a plasma 
([lo|). the increase of dissociation of weak electrolytes ( 11|) and the complex electrokinetics 
of charged colloidal spheres (Q|) under high electric flelds). It turns out that the motion 
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of a nanoparticle through a semidilute polymer is quite similar, though not identical, to 
that of convective diffusion occurring near a particle translating through a solution of some 
molecular species, which is continuously absorbed at the surface of the particle. Hence, 
Levich's analysis proves to be useful ([9]). The Peclet number Pe = va/D is important in 
the theory of convective diffusion. A sphere of radius a is surrounded by a zone of size a, 
partly depleted of absorbing, diffusing molecules. The time scale a^/D of diffusion where D 
is the diffusion coefficient of the molecules, has to be balanced against the time scale a/v 
of convection if the sphere has a speed v. At high Peclet numbers (Pe ^ 1), convection 
has a dramatic impact on the shape of the diffusional layer (j^). In the polymer case, 
the analogous phenomenon may be termed convective depletion but it is determined by a 
different dimensionless number defined by 

A = ^ (1) 

In the case at hand ([sll), the depletion layer of size a pertains to a section of polymer chain 
consisting of a'^/A'^ depleted segments of Kuhn Length A and diffusion coefficient Ds- I thus 
introduce eq (1) in the Rouse approximation. Accordingly, one expects the depletion layer 
to be deformed strongly by convection at high A. 



2 Convective depletion at high velocities 

The dynamics of a nanosphere in a semidilute solution of flexible polymer of concentration 

n 

Co has been discussed previously ([5]). The free energy of depletion F is solely of entropic 
origin and is a functional of c(T^,t), the density of polymer segments at position "r^ and 
time t 

F[c(r,t)] = ^A'kBT J dr (Vc^)' (2) 

Here, is Boltzmann's constant and T is the temperature. For convenience, the sphere is 
held fixed and its centre is the origin of our Cartesian coordinate system. At great distances 
from the sphere, the polymer solution flows uniformly with velocity —Vq in the z direction. 
The flow is perturbed by two effects: the solvent is assumed to stick to the sphere, and 
the polymer segments cannot penetrate its surface; the stick boundary condition seems 
reasonable since the polymer density tends to zero at the surface (i.e. we have depletion) 
and we may argue that the segments within the depletion layer are essentially Rouselike, 
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exerting little backflow on the solvent itself The radius a of the sphere is much 

smaller than both the polymer correlation and hydrodynamic screening lengths. Hence, the 
velocity ll (T^, t) of a test segment is given by a balance of forces 

Tt = Ij" + m f (3) 

where ^ {~f^,t) is the background velocity of the incompressible solvent, / is the force on 
the test segment by virtue of its connection to some chain section and m is the mobility of 
the segment. 

The gradient of the chemical potential V/i = — f now drives the reorganization of 
the segment distribution via the equation of continuity and eqs (j21) and 0; the chemical 
potential of a segment is fx {T^, t) = 5F/5c (^, t). 

^ + ^ • V {^^) = -BV ■ [v'V {^-'^v)) (4) 

I have here introduced B = ^niA^kBT and a dimensionless "wavefunction" Lp{~r^t) = 
^ 1 

{c{r ,t) /cq)^ . This equation has been solved for a sphere translating slowly under stationary 
conditions in the absence of convection (Q)- Here, I wish to investigate convective depletion 
for fast uniform propagation of the sphere in the steady state. It turns out to be expedient 
to let the Stokes stream function ip become the independent variable in problems involving 
convection {^) 

^ = -2^osm6'lr ~2^'"~^2~/ ^ ' 

-v,cose[l-- + — ] (6) 



sin 6 36 \ 2r 2r^ ^ 

ve = — ^ = -t;osme 2- — -— 7 

rsmO or 2 \ 2r 2r-^ J 

The Stokes stream function and velocity of the solvent, here given in terms of spherical 
coordinates (r, 6, ips) (Q]), are assumed to be unperturbed by the polymer within the zone 
of depletion. I next switch from the variables (r, 9) to the new independent variables (6', ip) 
assuming axisymmetry . It is readily shown that the convective term in eq (4) reduces to 



d r r 86 

with the help of eqs (6) and (7). 



(8) 



At this juncture, I introduce an approximation at high A. Levich has treated the con- 
vective diffusion of a fast moving sphere at high Peclet number (Q) and it is expedient to 
follow closely the spirit of his analysis of the thin boundary layer in front of such a particle. 
There is a stagnation point at = and r = 0, so at high velocities one may picture the 
region of polymer depletion as in fig. 1. It is compressed by convection into a layer thinner 
than the radius of the sphere. Hence, I introduce the coordinate y = r — a with y ^ a in 
the regime of interest and approximate eqs (5-7) to the leading order as follows 



(9) 



I d%l) 3 

— W - COSU 



(10) 



1 dip 3 fy\ . „ 

vg ~ ^-p:^ -w { -]smu 

a sm u oy 2 



(11) 



These equations describe the unperturbed flow of the incompressible solvent in the boundary 
layer. Eq (4) is now rewritten as 



Vg dif^ 



a 89 



-B 



d_ 

dy 



dy 



(12) 



with the help of eqs (8-11). Hence, we attain an expression for the convective depletion 
within the boundary layer at the front of the sphere as a function of two new independent 
variables: the Stokes stream function ip by first eliminating y using eqs (9) and (11), and 
the variable s for convenience. 



ds 



d_ 

dip 



-1 r^^^w 



dip W-i^- 



dip 



(13) 



s = Ba^ {'iwf'^l{e) + ci 



(14) 



I{9)= / dOsin^e (15) 

Here, ci is a constant. 
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It is now useful to introduce the boundary conditions. 



(f = at ijj — 



(16) 



ip — 1 a,s ijj 



— oo 



(17) 



(p = 1 at e = 0,ip = 



(18) 



&^(p d'^if dip 



at y = 



(19) 



e 



The first equation expresses the usual depletion of the polymer at the surface which I suppose 
to be independent of the velocity of the particle. Far away from it, the polymer density must 
tend to the bulk value (eq (17). However, eq (16) does not hold at ^ = for there is a point 
of stagnation at ^ = 0, = so we require the boundary condition given by eq (18). 
Finally, polymer segments cannot penetrate the sphere so the radial flux cu of polymer to 
the surface must be zero. This leads to the condition given by eq (19) upon temporarily 
reverting to the variable y. In view of the boundary conditions and the form of eq (13), it is 
now possible to posit a similarity solution (/? (77) in terms of the variable 77 = — The 



resulting equation is then further simplified with the help of the new independent variable 



This nonlinear equation looks rather awkward but a local analysis turns out to be useful. 
3 Analysis of eq (20) 

Usually one would substitute p = exp b and j = db/dp into eq (20) to decrease the order 
of the equation by one but this does not appear to be useful here. I have not attempted 
to derive a uniformly valid approximation to eq (20) because the term {(Pp/dp^)^ appears 
to play a significant role at p — O (1). Nevertheless, a local analysis within three regimes is 
fruitful for this term happens to be negligible close to and far away from the sphere. 




(20) 
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3.1 Inner region 



At the surface of the sphere, the polymer density vanishes (eq (16)) so that a series 
expansion of the form 

^ = C2 {p" + c^p"^" + c^p''^^^ + ..) (21) 

may be assumed with a > and /? > and where C2, C3 and C4 are constants. Upon 
substituting this into eq (20), one discerns that there are only two viable solutions, if a = 1 
and a = 3/2. Next if the radial polymer flux to the sphere is required to vanish (eq (19)), 
a must either be equal to unity or be greater than 3/2. Hence, the only feasible solution 
for is if a = 1. Again inserting eq (21) into eq (20) establishes the value of /3 = 5 and the 
constants C3 and C4 

The term in eq (20) containing the second derivative contributes only at order p^^ in eq (22). 
Next, the boundary condition at the point of stagnation (eq 18) determines the value of the 
constant ci = 0; at ^ = 0, p must be nonzero for y > 0. 

3.2 Outer region 

Next, I attempt a WKB approximation in the region far away from the sphere by writing 
(p = 1 + h with \h\ <^ l.The term containing the second derivative is then subdominant and 
the order of the resulting linear equation is reduced by unity with the aid of the substitution 

g = dh/dp 

There seems to be a turning point at p = but I focus on the outer region p ^ 1 only, 
since there is at least one intermediate regime centered around p = O (JQ where {d'^cp/dp'^)'^ 
competes with ipd^ip/ dp^ . There are three WKB solutions to eq (23) ( 14l |) 

' exp So (24) 



S, = wldp[^^] (25) 
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where 1, — | + |-\/3i and — | — ^VSi are the three roots of the equation — l.The outer 
solution obeying eq (17) is then 

(/p^l + cs^ dp p-h-'^^P'^' sin (^^VSHp^/^ + ce^ (26) 

5 V 5 

3.3 Intermediate region 

According to cqs (22) and (26), the (square root of) the polymer density exhibits the 
oscillatory profile depicted in fig (2) with an as yet unknown intermediate regime at around 
p — O {1) with at least one maximum, the first, centered at pi say. We Taylor expand (p (p) 
around pi 

(fip) ^(f (pi) - ^cr (p - pif + ^cg {p - pif + ^cg {p - Pi)'' (27) 

Upon substituting this in eq (20), we see that the coefficient cg is positive as it should be. 
One expects pi = O (1) and (fipi) — O {1) so that cy and Cg are also 0(1) for consistency. 

It is difficult to ascertain the stabihty of the solution determined by eq (20) with regard 
to eq (4) under general circumstances in view of the inhomogeneity of the fiuid velocity 
coupled to the nonlinear depletion (even though it is quasilinear) . Nevertheless, it can be 
shown that the solution is Hadamard stable. 



I/O 



4 Force on the sphere 

Let us first rewrite the similarity variable p in terms of the original coordinates 9 and y 
pertaining to the boundary layer 

1 f 3v \ usin^ 

Accordingly, I introduce an angular-dependent thickness of the layer defined by 

\ V J sm.9 

Both p and 5 (9) are well behaved as 9 tends to zero and valid right up to ^ = 7r/2, to a 
good approximation. I note that they are independent of the bulk density Cq because eq (4) 
is quasilinear. 



Next, the force on the sphere is defined by 

1 



/ 



dr c{r ) 



V/i 



-A coksT I dr if 



(30) 



The integration apphes only to the region definable as the boundary layer. Since 5q a,l 
introduce a Derjaguin approximation where the distance of a point with coordinates (i?, y) to 



the surface of the sphere is given by y — a + a/ {a? + E?) ~ ?/ + (i?^/2a).The radial coordinate 
R is perpendicular to the z axis. The integral in eq (30) is thus rewritten as 



pa pa 

J~ / dR2nR dy G [y + [R^ /2a)) 
Jo Jo 

d^(p dip d'^cp 



(31) 



G (y) ^ V- 



0) 



dy'^ dy dy"^ 

To the leading order, J may be simplified after several integrations by parts and the 
substitution q = R^/2a 
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poo poo pc 

J ~ 27ra / dq dy G {q + y) = Ana 

Jo Jo Jo 



dy 



dip 
dy 



1 61=0 



6{0) 



(32) 



We conclude that the force on the sphere is inversely proportional to the thickness of the 
boundary layer 6 (0). 



/ 



A^acokBT/5 (0) 



(33) 



This is consistent with the force of depletion 



/ 



a^CQv/m = A'^coAkBT at low velocities 



derived previously ((^]); the crossover is at A = O (1) and 6 (0) ~ a. 



5 Concluding remarks 

The force on a sphere arising from convective depletion is illustrated for a globular protein 
propagating fast through a semidilute polymer under the influence of a constant uniform 
electric field E. If the protein bears Z positive elementary charges of charge g, the force 
ZqE on the protein ultimately causes it to translate uniformly with a velocity "i? given by 
eqs (29 ) and (33), 



DsA^ 


( ZqE \ 




f zn \ 


a3 


[A^cokBTj 







(34) 



Here it is convenient to introduce VL = qEQ/ksT, a dimensionless electric field scaled by 
the Bjerrum length Q{Q = 0.71 nm for water at room temperature; Q = 2.76 x lO^^E if 
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E is given in units of V / cm). Eq (34) is only valid at high enough fields: E > with 
^ A^Qcq/Z because the crossover must occur when the dimensionless quantity given 
by eq (1) is of order unity; E^, is about 1 kV / cm for a protein bearing 10 positive charges 
translating through a polymer solution of volume fraction 0.1. The electrophoretic velocity 
is thus a sensitive function of the protein charge and its dimensions. Electrophoresis in 
polymer gels at high electric fields could well be a useful technique for the separation of a 
mixture of proteins. Very high fields are currently employed in microchip electrophoresis 
because Joule heating may be kept to a minimum ([isj]). 

Figure Captions 

Fig. 1 Fixed sphere in a fiuid moving uniformly with velocity — fo in the z direction at 
infinity. The dashed curve demarcates approximately the onset of the depletion layer. 
Fig. 2 Square onset of the scaled polymer density as a function of the similarity variable 

p. 
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